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Q I Recent experiments on production of AA Hypernuclei have stimulated renewed interest 

■ in extracting the AA interaction from the few events identified since the inception of this 
field forty years ago. Few-body calculations relating to this issue are reviewed, particularly 
with respect to the possibility that A = 4 marks the onset of AA binding to nuclei. The 
Nijmegen soft-core model potentials NSC97 qualitatively agree with the strength of the 

^ I AA interaction deduced from the newly determined binding energy of Applying the 

O I extended NSC97 model to stranger nuclear systems suggests that A = 6 marks the onset 
^ ' of S binding, with a particle stable aJHc, and that strange hadronic matter is robustly 

■ bound. 
O _ 

1. INTRODUCTION 

O I Until 2001 only three candidates existed for AA hypernuclei observed in emulsion exper- 
^ I iments [1-3]. The AA binding energies deduced from these emulsion events indicated that 
the AA interaction is strongly attractive in the ^5*0 channel [4-6], with a AA pairing energy 
. - ASaa ~ 4.5 MeV, although it had been realized [7,8] that the binding energies of \^Be [1] 
! and AA-^e [2] are inconsistent with each other. This outlook has undergone an important 
\ change following the very recent report by the KEK hybrid-emulsion experiment E373 
of a well-estabhshed new candidate [9] for AAHe, with binding energy (Ai?AA ~ 1 MeV) 
substantially lower than that deduced from the older, dubious event [2]. Furthermore, 
there are also indications from the AGS experiment E906 for the production of light AA 
hypernuclei [10], perhaps as light even as aaH, in the {K~,K~^) reaction on ^Be. 

Since data on hyperon-nucleon [YN) and hyperon-hyperon {YY) interactions are scarce 
or even not readily available from laboratory experiments, the study of multistrange sys- 
tems provides a fairly exclusive test of microscopic models for the baryon-baryon (BB) 
interaction. The Nijmegen group has constructed over the years a number of one-boson- 
exchange (QBE) models (reviewed by Rijken in Ref. [11] and in these proceedings) for the 
BB interaction using SU(3)-flavor symmetry to relate coupling constants and phenomeno- 
logical short-distance hard or soft cores. In all of these rather different BB interaction 
models only 35 YN low-energy, generally imprecise data points serve the purpose of steer- 
ing phenomenologically the extrapolation from the A^A^ sector, which relies on thousands 
of data points, into the strange YN and YY sectors. It is therefore of utmost impor- 
tance to confront these models with the new AA hypernuclear data in order to provide 
meaningful constraints on the extrapolation to strangeness S = —2 and beyond. 
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Figure 1. Nijmegen OBE phase-equivalent 
soft-core AA potentials [12]. 



Figure 2. s-wave Faddeev calculations [12] 
of ASAA(A;^He) vs. A5AA(AiH, ^^He). 



The YN and YY s-wave interaction input potentials to the structure calculations here 
reviewed often consist of combinations of Gaussians with different ranges, such as to make 
these single-channel potentials phase equivalent to the Nijmegen OBE-model coupled- 
channel potentials. Of the several A A potentials due to Nijmegen models which are 
shown in Fig.l, NSC97c is the weakest one, of the order of magnitude required to reproduce 
A^aaIaaHc). The AA interaction is fairly weak for all six versions (a)-(/) of the Nijmegen 
soft-core model NSC97 [13], and versions e and / provide a reasonable description of 
single-A hypernuclei [14]. 



2. AA HYPERNUCLEI 

In this section I will review topical theoretical work on some of the light AA hypernuclear 
species connected to old and to new experiments. The anticipated existence of AA^e, now 
solidly established also experimentally [9], leads one to enquire where the onset of AA 
binding occurs. It was argued long ago that the three-body AAA^ system is unbound [15], 
and hence I will concentrate on the A — AA hypernuclear systems. Among the few 
heavier species reported todate, AA^e will be discussed briefly. 



2.1. ^/H - ^^He 

Figure 2 demonstrates a nearly hnear correlation between Faddeev-calculated values of 
ASaaIaaHg) and ABaaIaaH, aaHg), using several AA interactions which include (the 



lowest-left point) Vaa — [12]. Here 
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Figure 3. s-wave FY calculations [16] for 
Apn, AAd and AApn. 



Figure 4. s-wave Faddeev calculations of 
^aa(aaH) in a AAd model [16]. 



where B/^/^{j^Z) is the AA binding energy of the hypernucleus j^Z and Bj^('\^^Z) is 
the (2J+l)-average of B\ values for the ^"^^^^ Z hypernuclear core levels. ASaa increases 
monotonically with the strength of Vaa, starting in approximately zero as Vaa 0, 
which is a general feature of thrcc-body models such as the aAA, ^HAA and ^HeAA 
models used in these s-wave Faddeev calculations [12], and also as shown below for dAA 
s-wave Faddeev calculations [16]. The 7=1/2 ^^H - AA^e hypernuclei are then found 
to be particle stable for all the AA attractive potentials here used. This conclusion holds 
also when the s-wave approximation is relaxed [17]. 



2.2. ^/H 

I start by discussing the first Faddeev- Yakubovsky (FY) four-body calculation of ^^^H 
[16]. For two identical hyperons and two essentially identical nucleons (upon introducing 

isospin) as appropriate to a AApn model calculation of ^aH, the 18 FY components 
reduce to seven independent components satisfying coupled equations. Six rearrangement 
channels are involved in the s-wave calculation [16] for ^^H(l+): 



{ANN)s=r+A, {ANN)s=s+A, {AAN)s=.+N 
for 3-1-1 breakup clusters, and 

(AA)5=o + {NN)s=i , {AN)s + (AiV)^, 



(2) 



(3) 



with {S, S') = {0, 1) + (1, 0) and (1, 1) for the latter 2+2 breakup clusters. 

Using Vaa which reproduces 5AA(AAHe), the four-body calculation converges well as 
function of the number N of the FY basis functions allowed in, yet it yields no bound 
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Figure 5. Saa(aaH) calculated by Ne- 
mura et al. [19] using the stochastic vari- 
ational method. 



Figure 6. Saa(aaH) calculated in a kkd 
model, using exponential Ac? potentials 
constrained by NSC97f and NSC97f' [16]. 



state for the AApn system, as demonstrated in Fig.3 by the location of the 'AApn' curve 
above the horizontal straight line marking the 'A + threshold'.^ In fact these FY 
calculations exhibit little sensitivity to Vaa over a wide range. Even for considerably 
stronger AA interactions one gets a bound ^^H only if the AA^ interaction is made con- 
siderably stronger, by as much as 40%. With four AA^ pairwise interactions out of a total 
of six, the strength of the AA^ interaction (about half of that for A^A^) plays a major 
role in the four-body AApn problem. However, fitting a Ad potential to the low-energy 
parameters of the s-wave Faddeev calculation for Apn, and solving the s-wave Faddeev 
equations for a AAd model of aa^, this latter four-body system is calculated to yield a 
1+ bound state, as shown in the figure by the location of the asymptote of the 'AAd' 
curve below the 'A -|- threshold'. The onset of particle stabihty for ^^H(l+) requires 
then a minimum strength for Vaa which is exceeded by the choice of SAA(AAHe) [9] as 
a normalizing datum (equivalent to — oaa ~ 0.8 fm [12]). This is demonstrated in Fig.4 
where Faddcev-calculated Bj^^^j^^B.) values are shown as function of the AA scattering 
length oaa for two different functional forms of the fitted Ad potential. Disregarding spin 
it can be shown that, for essentially attractive AA interactions and for a static nuclear 
core d, a two-body Ad bound state implies binding for the three-body AAd system [18]. 
However, for a nan static nuclear core d (made out of dynamically interacting proton and 
neutron), a Ac? bound state does not necessarily imply binding for the AAc? system. It is 



^This threshold was obtained as the asymptote of the Apn s-wavc Faddeev calculation which uses model 
NSC97f [14] for the underlying AN interaction, yielding Ba{IR{^^)) = 0.19 MeV. Using model NSC97e, 
with B\{\}i{^^)) = 0.07 MeV, does not alter the conclusions listed below. 
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Figure 7. s-wave FY calculations [12] for Figure 8. NSC97 phase-equivalent YY po- 
/^Be: ^BeAAvs. aaAA. tentials [20]. 



questionable whether by incorporating higher partial waves, or AA — 5A^ coupling effects, 
this qualitative feature will change. 

The above conclusions have been very recently challenged by Nemura et al.[19]. Fig. 5 
demonstrates that within their stochastic variational calculation, which uses the NSC97f 
input of the Filikhin and Gal calculation [16] for the various pairwise interactions, a 
'pnAA' model always yields more binding than a 'rfAA' model does. Particle stability 
for ^^H(l"'") in this variational calculation requires a minimum strength for Vaa which 
is exceeded by the choice of B\a{aa^q) [9] as a normalizing datum. Yet, Nemura et al. 
argue that the AN interaction in the channel, when adjusted to the binding energy 
calculated for the A = 4 A hypernuclei, should be taken weaker than that used by Filikhin 
and Gal and that, when this constraint is implemented ('set A' in the figure), particle 
stability for ^^H(l+) requires a minimum strength for V\\ which is not satisfied by the 
choice of Bw^j^^Re) as a normalizing datum. A similar strong dependence on the ^^S" 
AA^ interaction within a AAd model is shown in Fig. 6 for two versions / and /' of model 
NSC97 which produce the same B\{\}i{^^)), while differing slightly by the location of 
the spin-flip |^ excited state [16]. 



2.3. 

For heavier AA hypernuclei, the relationship between the three-body and four-body 

models is opposite to that found by Filikhin and Gal for aaH: the AAC1C2 calculation 
provides /^i^f/ier binding than a properly defined AAC calculation yields (with C = C1 + C2) 
due to the attraction induced by the AC1-AC2, AAC1-C2, C1-AAC2 four-body rearrange- 
ment channels that include bound states for which there is no room in the three-body 
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Figure 9. s-wave Faddeev calculations [20] for j^^B. and aIHc. 



AAC model. The binding energy calculated within the four-body model increases then 
'normally' with the strength of Vaa [12]. This is demonstrated in Fig. 7 for aaBc using 
several AA interactions, including Vaa = which corresponds to the lowest point on each 
one of the straight lines. The origin of the dashed axes corresponds to Ai?AA = 0. Within 
the 4-body aaAA model, the fairly large value AB\\{j^^Be) ~ 1.5 MeV in the limit 
Vaa ^ is due to the special aa cluster structure of the ^Be core. The correlation noted 
in the figure between a^a^c and aa^c calculations, and the consistency between various 
reports on their Baa values, are discussed by Filikhin and Gal [12,20]. In particular, the 
two solid points next to the lowest one on the '4-body model' line in Fig. 7, corresponding 
to two versions of model NSC97 [13], are close to reproducing (the 'new') SaaIaaHc) but 
are short of reproducing (the 'old') Baa{aa^G:) by about 2.3 ±0.4. This apparent discrep- 
ancy may be substantially reduced by accepting a a^a^c weak decay scheme that involves 
the 3 MeV excited aBc doublet rather than the aBc ground state [21]. This conclusion 
may also be inferred from the recent 4-body calculations by Hiyama et al. for A = 7 —10 
AA hypernuclei [22]. 

3. THE ONSET OF S STABILITY 

Since model NSC97 [14] provides a qualitatively successful extrapolation from fits to 
A'^A^ and YN data to 5" = — 2, and noting the strongly attractive ^5*0 AS potentials shown 
in Fig. 8 in comparison to the fairly weak AA potentials when model NSC97 is extrap- 
olated to S = —3, —4 [13], it is natural to search for stability oi A = Q, S = —3 systems 
obtained from aaHc upon replacing one of the A's by S. Faddeev calculations [20] for the 
0"*" / = 1/2 ground-state of a|H and aIHc, considered as q;A5~ and q;AS° three-body 
systems respectively, indicate that aIHc is particle-stable against A emission to aa^c for 
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potentials simulating model NSC97, particularly versions e and /, whereas is unsta- 
ble since M(S~) > M(S°) by 6.5 MeV.^ This is demonstrated in Fig. 9. Nevertheless, 
predicting particle stability for ^|He is not independent of the assumptions made on the 
experimentally unexplored Ea interaction which was extrapolated from recent data on 
^^C [23]; hence this prediction cannot be considered conclusive. 

4. STRANGE HADRONIC MATTER 

Bodmer [24], and more specifically Witten [25], suggested that strange quark matter, 
with roughly equal composition of u, d and s quarks, might provide the absolutely stable 
form of matter. Metastable strange quark matter has been studied by Chin and Kerman 
[26]. Jaffe and collaborators [27,28] subsequently charted the various scenarios possible 
for the stability of strange quark matter, from absolute stability down to metastability 
due to weak decays. Finite strange quark systems, so called strangelets, have also been 
considered [27,29]. 

Less advertised, perhaps, is the observation made by Schaffner et al. that metastable 
strange systems with similar properties, i.e. a strangeness fraction fs = —S/A k, 1 and 
a charge fraction fq — Z/A might also exist in the hadronic basis at moderate 

values of density, between twice and three times nuclear matter density [30,31]. These 
strange systems are made out of iV, A and S baryons. The metastability of these strange 
hadronic systems was established by extending relativistic mean field (RMF) calculations 
from ordinary nuclei {fs = 0) to multi-strange nuclei with fs 0. Although the detailed 
pattern of metastability, as well as the actual values of the binding energy, depend specif- 
ically on the partly unknown hyperon potentials assumed in dense matter, the predicted 
phenomenon of metastability turned out to be robust in these calculations [32]. 

Recently, model NSC97 and its extension [13] were used to calculate within the RMF 
framework the minimum-energy equilibrium composition of bulk strange hadronic matter 
(SHM) made out of the SU(3) octet baryons N, A, S and S, over the entire range of 
strangeness fraction < fs < 2 [33]. The main result is that SHM is comfortably 
metastable in this model for any allowed value of fs > 0. The iVAS composition and the 
binding energy calculated for equilibrium configurations with < 1 resemble those of 
model 2 in Refs. [30,31]. The extension of model NSC97 [13] yields particularly attractive 
SS, EE and ES interactions, but fairly weak AA and NS interactions. Consequently, for 
/5 > 1, S's replace A's due to their exceptionally strong attraction to S and S hyperons. 
As is shown below, a first-order phase transition occurs from A^AS dominated matter for 
/5 < 1 to A^ES dominated matter for > 1, with binding energies per baryon reaching 
as much as 80 McV. 

A phase transition is visualized in Fig. 10 where the binding energy is drawn versus the 
baryon density for several representative fixed values of fs- For fs — 0.8, there is a global 
minimum at a baryon density of pb = 0.27 fm~^. A shallow local minimum is seen at 
larger baryon density at pB = 0.72 fm~^. Increasing the strangeness fraction to fs = 0.9 
lowers substantially the local minimum by about 20 McV, whereas the global minimum 
barely changes. At fs = 1.0 this trend is amplified and the relationship between the 

^Recall that the 7=1/2 ^^H - aaHb hypernuclei, within a AAC Faddeev calculation, are particle stable 
even in the limit Vaa 0. 
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Figure 10. Transition from A^AS to NEE 
matter upon increasing the strangeness 
fraction [33]. 



Figure 11. Strange hadronic matter com- 
position as function of the strangeness 
fraction [33]. 



two minima is reversed, as the minimum at higher baryon density becomes energetically 
favored. The system will then undergo a transition from the low density state to the 
high density state. Due to the barrier between the two minima, it is a first-order phase 
transition from one minimum to the other. 

Fig. 11 demonstrates explicitly that the phase transition involves transformation from 
A'"AS dominated matter to A^SS dominated matter, by showing the calculated composi- 
tion of SHM for this model (denoted N) as function of the strangeness fraction fs- The 
particle fractions for each baryon species change as function of fs- At = 0, one has 
pure nuclear matter, whereas at fs — 2 one has pure S matter. In between, matter is 
composed of baryons as dictated by chemical equilibrium. A change in the particle frac- 
tion may occur quite drastically when new particles appear, or existing ones disappear in 
the medium. A sudden change in the composition is seen in Fig. 11 for fs = 0.2 when 
S's emerge in the medium, or at fs = 1.45 when nucleons disappear. The situation at 
fs — 0.95 is a special one, as E's appear in the medium, marking the first-order phase 
transition observed in the previous figure. The baryon composition alters completely at 
that point, from NE baryons plus a rapidly vanishing fraction of A's into EH hyperons 
plus a decreasing fraction of nucleons. At the very deep minimum of the binding energy 
curve (Fig. 3 of Ref. [33]) SHM is composed mainly of E's and S's with a very small 
admixture of nucleons. 
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5. CONCLUSION 

I have presented Faddeev calculations for j^^R - ^^He and j^^He, and first ever four- 
body Faddeev-Yakubovsky calculations for ^yaH and aaBc, using two-body potentials 
fitted to the low-energy scattering parameters or to the binding energies of the respective 
subsystems. In particular, for ^^H, A^A^ and AA^ interaction potentials that fit the binding 
energy of were used. No ^^^H bound state was obtained for a wide range of AA 
interactions, including that corresponding to SAA(AAHe). This non binding is due to the 
relatively weak AA^ interaction, in stark contrast to the results of a 'reasonable' three- 
body AAd Faddeev calculation. Further experimental work is needed to decide whether or 
not the events reported in the AGS experiment E906 [10] correspond to ^aH, particularly 
in view of subsequent conflicting theoretical analyses [34,35]. More theoretical work, 
particularly on the effects of including explicitly AA— SA?"— EE channel couplings, is called 
for. Preliminary estimates for such effects within the NSC97 model, or its simulation, 
have been recently made [17,36-39]. In particular, Lanskoy and Yamamoto [39] have 
focused attention to the substantial charge symmetry breaking effects introduced by the 
different S— ^He thresholds into the binding energy calculation of the ^^H - j^lUe ground 
states. In addition to increasing the calculated ASaa value of aaHg with respect to that 
of aaH, on top of the difference already shown in Fig. 2, it is found that for model 
NSC97 Ai?AA(AAHc) gets as large and perhaps even larger than that of ^vaHc (~ 1 MeV). 
A simultaneous determination of the binding energies of AA^e and aaHc would help to 
discriminate between several versions of QBE models which differ markedly from each 
other regarding the strength of the off-diagonal AA — EN couphng. 

Accepting the predictive power of model NSC97, Faddeev calculations suggest that 
^iHe may be the lightest particle-stable S = —3 hypernucleus, and the lightest and least 
strange particle-stable hypernucleus in which a S hyperon is bound. Unfortunately, the 
direct production of AH hypernuclei is beyond present experimental capabilities, requiring 
the use of Q~ initiated reactions. 

Finally, I have focused on the consequences of using model NSC97 for the binding and 
composition of strange hadronic matter. Strange hadronic matter is comfortably stable, 
up to weak decays, over a wide range of baryon-baryon interaction models, including 
model NSC97 here chosen because it successfully extrapolates from the 5" = 0, — 1 sectors 
in which it was constructed into the S — —2 sector, nearly reproducing EaaCaaHc). The 
phase transition considered in this review has been recently discussed by the Frankfurt 
group [40] in the context of phase transition to hyperon matter in neutron stars. Un- 
fortunately, it will take lots of imagination to devise experimentally a way to determine 
how attractive those AS, SS, SE, EE interactions are, which are so crucial for the results 
exhibited in this review. 
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